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ABSTRACT 

We demonstrate how the observed H2O ice column densities toward three dense glob- 
ules in the Southern Coalsack could be used to constrain the ages of these sources. 
We derive ages of ~ 10 5 yr, in agreement with dynamical studies of these objects. We 
have modelled the chemical evolution of the globules, and show how the molecular 
abundances are controlled by both the gas density and the initial chemical conditions 
as the globules formed. Based on our derived ages, we predict the column densities 
of several species of interest. These predictions should be straightforward to test by 
performing molecular line observations. 

Key words: astrochemistry - ISM: clouds - ISM: individual objects: Coalsack - 
ISM: molecules - molecular processes 



> 

cn 

o 
o 



x 



1 INTRODUCTION 

Chemical reactions alter the composition of dense interstel- 
lar gas and theoretical models of this evolution offer the 
prospect of determining the relative ages of dense cores and 
molecular clouds (e.g., Stahler 1984; Leung, Herbst & Hueb- 
ner 1984). The accuracy of these derived ages will increase 
with that of the initial conditions of the model calculations, 
i.e. knowledge of the chemical state of the gas at some ar- 
bitrary time in the past. In regions of active star formation 
the prompt (£ <, lyr) removal of icy grain mantles of known 
composition, either by thermal evaporation or sputtering in 
shocks, almost simultaneously alters the chemical compo- 
sition of protostellar cores and leads to a prolonged phase 
(~ 10 years) of 'hot core' chemical evolution. The destruc- 
tion of abundant 'parent' molecules and the formation of 
'daughter' species in principle allows the time since man- 
tle removal to be measured from the presently observed gas 
phase composition (Charnley, Tielens & Millar 1992; Charn- 
ley 1997; Hatchell et al. 1998). 

Despite the fact that they are probably formed from 
more diffuse material, of fairly simple chemical composition, 
estimating relative 'chemical ages' for molecular clouds and 
for dense (starless) cores within those clouds is more prob- 
lematic. Models of quiescent dark clouds have largely been 
based on so-called 'pseudo-time-dependent' models (e.g., Le- 
ung et al. 1984). In this case, the physical conditions re- 
main fixed at dense cloud values throughout the calculation 
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whereas the initial conditions for the heavy elements are 
those expected from the diffuse interstellar medium. The 
chemical evolution in these purely gas phase models can be 
roughly divided between the abundances present at 'early- 
times', prior to a major fraction of the gas phase carbon 
being incorporated into CO (~ 10 5 years), and 'late-times', 
approximately when the chemistry has reached a steady- 
state (~ 10 6 - 10 7 years). Time-dependent models of dark 
clouds have shown that the carbon-chain molecules are best 
explained by 'early-time' chemistry whereas several other 
molecular abundances, such as N2H + and NH3, are closer to 
the predicted steady-state values. This has led to the idea 
that observed spatial compositional gradients (e.g. involv- 
ing the carbon-chains and ammonia) are due to a differing 
chemical ages in individual clouds (e.g. Hirahara et al. 1992; 
Taylor, Morata & Williams 1998). 

Much effort has been made by theoretical astrochemists 
to reconcile the apparent chemical youth of molecular clouds 
with other evidence which suggests that these objects are 
long-lived entities (e.g., Nejad, Williams & Charnley 1990). 
For example, in the 'classical' picture of molecular cloud evo- 
lution and low-mass star formation (Shu, Adams & Lizano 
1987), the ambipolar diffusion time-scale of a plasma with 
fractional ionization x c 
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(McKee et al. 1993) governs the loss of magnetic support 
which is a prerequisite for gravitational collapse (Shu et al. 
1987). In this picture, molecular clouds are inferred to have 
lifetimes of ~ 10 Myr (Mouschovias, Tassis & Kunz 2006). 



© 0000 RAS 



2 S.D. Rodgers et al. 



However, recent studies suggest that molecular cloud forma- 
tion, dissipation and star formation, are all rapid dynamical 
processes (Elmegreen 2000; Hartmann, Ballasteros-Paredes 
& Bergin 2001; Mac Low & Klessen 2004). In this revised 
picture, clouds form through the dissipation of supersonic 
turbulence and the important dynamical time-scale is the 
gravitational free-fall time at number density riH 

tg = 4.35 x 10 7 rc H 1/2 yr (2) 

Although the debate over the mode of star formation is 
not yet settled (e.g., Mouschovias et al. 2006; Ballesteros- 
Paredes & Hartmann 2006), if dense clouds and stars form 
within a few free-fall times, with the cloud life-time being ~ 
1-3 Myr (Hartmann et al. 2001), determining the chemical 
ages of molecular clouds may, after all, only require consid- 
eration of the abundances during the 'early time' phase. 

To test and refine the idea that chemical models could 
be used to measure the (relative) ages of molecular clouds re- 
quires two additional criteria. First, we would like to model 
the chemical evolution in a molecular cloud for which there 
is independent evidence that it formed recently. The South- 
ern Coalsack is a well-studied region which fulfills this cri- 
terion. The Coalsack is a group of molecular globules (or 
clumps) where the low average gas densities and the lack 
of discernible star formation activity (Nyman, Bronfman & 
Thaddeus 1989; Kato et al. 1999), the presence of a dif- 
fuse inter-clump medium, as evidenced by CH + observations 
(Centurion, Cassola & Vladilo 1995), and the apparent dy- 
namical youth of the major globule (Lada et al. 2004), all 
suggest that it is a region where simplistic chemical models 
may actually be applicable. Second, we would like to use 
another time-scale to constrain the chemical evolution. The 
time-scale over which icy grain mantles are laid down in 
molecular clouds provides a natural constraint that is inde- 
pendent of the precise chemical state of the gas (Charnley, 
Whittet & Williams 1990). This time-scale, ti cc , is deter- 
mined by the rate at which oxygen atoms collide with and 
stick to dust grains, and depends on the gas density, the 
atomic oxygen abundance, and the total effective grain sur- 
face area (see Sj3j. 

Ice mantles have been detected towards several glob- 
ules in the Southern Coalsack by Smith et al. (2002) and 
this allows us to estimate the mantle formation time-scale 
for each. In this paper, we use the mantle formation time- 
scale as the basic clock by which to measure the chemical 
chronology of dense globules in the Southern Coalsack, and 
thereby predict their relative chemical compositions. 



2 THE SOUTHERN COALSACK 

The Southern Coalsack is a conspicuous dark molecular 
cloud complex in the southern sky. Unlike similar regions 
in the northern Milky Way, such as Taurus, Ophiuchius 
and Perseus, a striking feature of the Southern Coalsack 
is the, apparently complete, absence of star formation (Ny- 
man et al. 1989; Kato et al. 1999; Bourke, Hyland & Robin- 
son 1995; Racca, Gomez & Kenyon 2002). Optical and 
millimetre- wave observations of molecules towards the back- 
ground field star HD 110432, which sample a line-of-sight 
through the inter-clump medium, show that the Southern 
Coalsack is pervaded by translucent (Ay ~ 1 — 3 mag.), 



diffuse (riH ~ 10 2 cm 3 ) material (Codina et al. 1984; van 
Dishoeck & Black 1989; Gredel, van Dishoeck & Black 1994; 
Centurion et al. 1995; Rachford et al. 2001). The low CO 
abundance (Codina et al. 1984), and the presence of only 
simple molecules, are consistent with the chemistry expected 
in translucent clouds. 

Tapia (1973) identified three dark globules (hereafter 
Globules 1,2 & 3) which showed higher than average extinc- 
tion. Apart from isotopomers of CO, the Southern Coal- 
sack has thus far only been mapped in OH, H2CO, and H I 
(Brooks, Sinclair & Manefield 1976; McClure-Grifnths et al. 
2001). Nyman et al. (1989) mapped the Southern Coalsack 
in CO lines and showed that Tapia's globules were also peaks 
in CO emission. Observations of 13 CO and C ls O have identi- 
fied several other distinct cores (Vilas-Boas, Myers & Fuller 
1994; Kato et al. 1999). These studies indicate that dense 
condensations in the Coalsack have densities in the range 
nu ~ 10 3 — 10 4 cm -3 with gas and dust temperatures of 
around 10 K (Jones, Hyland & Bailey 1984). An infrared 
extinction study of Globule 2 (Lada et al. 2004) indicates 
that it is in a state of dynamical equilibrium, and suggest 
that it formed only recently. Thus, the Southern Coalsack 
is a region where dense cloud cores have recently condensed 
out of a more diffuse, low-extinction, interclump medium. 

Smith et al. (2002) performed 3^tm spectroscopy to- 
wards eight field stars, located behind the Southern Coal- 
sack, to search for H2O ice absorption towards several of 
the dense globules. Six of these sources lay behind or near 
Globule 2 and one behind each of the other two. Smith et al. 
detected water ice in each of Tapia's globules with an addi- 
tional four tentative detections in Globule 2. The largest ice 
column densities were found toward stars with the largest 
visual extinctions, as found in many other interstellar clouds 
(e.g., Whittett et al. 1988). In general, the fact that the wa- 
ter features are only apparent above a certain level of Ay 
is not direct proof that the ice must be present in the dens- 
est regions, although it is suggestive that this is the case. 
In the Coalsack, however, the fact that Smith et al.'s field 
stars lay behind the globules, and we know from the CO 
line maps that these are regions of higher density, leads us 
to conclude that the observed H2O ice is frozen out in the 
cold, dense globules, and not in the inter-clump medium. In 
what follows, we use data from Smith et al. for their sources 
designated SS1-2, SS2-2 and SS3-1 to derive a mantle for- 
mation time-scale in each of Tapia's Globules 1, 2 & 3. 



3 ICE FORMATION TIME-SCALES 

We now show how observations of the ice column density can 
be used to estimate the ice mantle formation time-scales. 
Accretion and hydrogenation of oxygen atoms will lead to a 
spherical (silicate) grain core of radius a being covered by a 
water ice mantle, and the mean particle radius growing as 

r{a,t) — a + d(a,t) (3) 

where d(a, t) is the instantaneous mean mantle thickness. A 
spherical grain of instantaneous radius r(a,t) increases in 
size by accretion of mass through oxygen atoms sticking to 
the surface at a rate (Wickramasinghe 1967) 

dr S fkTmo\ 1/2 m fA\ 

Tt = 7- {-*r) no{t) (4) 
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where T is the gas temperature, 5* is the sticking efficiency, 
mo is the mass of an oxygen atom, pi cc is the buik density 
of water ice, and no(t) is the number density of O atoms in 
the gas. For simplicity, we have ignored the increase in mass 
due to the accretion of two hydrogen atoms per oxygen. Be- 
cause both the molecular accretion rate and the number of 
molecules required to build up one monolayer of coverage 
are proportional to the surface area of the grain, the ra- 
dial growth of grains is independent of grain size. Hence, at 
any particular time all grains will be covered by the same 
thickness of ice. 

An ice mantle will contain on average W(a,t) water 
molecules 



W(a,t) = 



54mH 



\a,t) 



(5) 



where ran is the hydrogen nucleon mass. A pathlength L 
through a molecular cloud with uniform dust and gas num- 
ber densities, n gr (a) and hh, has an ice column density of 



N ice (t) 



W(a, i)n gr (a) da dl 



(6) 



If the hydrogen column density, A?h (= nuL), is related to 
the visual extinction through 

(7) 



Nn = 1.6 x W 21 A V cm" 2 



(Bohlin, Savage & Drake 1978) then the time to accumulate 
this much ice, the mantle formation time-scale, tice, satisfies 



r 3/ , \ 3i -3.5 , p o w , n -20 Ni ce (ti ce 

[r (a, Uce) — a \a da — 5.8 x 10 — 



(8) 



Here we have assumed that the dust obeys an MRN size 
distribution (Mathis, Rumpl & Nordsieck 1977) such that 



r (a) = Aa 



n H 



(9) 



in a range bounded by a_ = 0.005 /im and a+ = 0.25 /im, 
and where A = 7.76 x 10~ 26 cm 2 ' 5 /H for silicate dust (Draine 
& Lee 1984). The ice mantle thickness at this time, d m = 
rf(iice), is just r(o,ti ce ) — a, and so we can evaulate the in- 
tegral in © to obtain a cubic equation in d m . As Nice/ Ay 
is determined observationally, the cubic can be solved to 
obtain d m for any source. 

We must now calculate the time, tice, required to accrete 
this thickness of ice. Accretion on to grains removes oxygen 
atoms from the gas at a rate 

z-^-rmS no(t) (10) 

4 \-n-m J 



dno 
dt 



where the total available grain surface area per nucleon at 
time t, Ag r (t), is 

A g r(i) = (47rr 2 (a, t)n gI (a)/ n H ) 



4nr (a,t)n gI (a)/nu da 



9.75 x 10" 



r (a,t)a 



da 



(11) 



The factor of 4 in eqn. (|10[) is due to the fact that A gr repre- 
sents the total grain surface area (47rr 2 ), whereas the cross- 
section for accretion onto each grain is equal to -irr 2 . 

Subject to some approximations, an analytical solution 
for ticc can be obtained. Firstly, accretion on to grains is 



the dominant loss route for O atoms in dark clouds, so we 
neglect the effects of other chemical processes on no(t). We 
also assume that the underlying power-law for n gr (a) is not 
modified by coagulation. For time-scales sufficiently short 
that Agr(t) is approximately unchanged (i.e. only thin man- 
tles have formed), the total surface area of the bare grains 
integrated over the MRN size distribution, A gr (0), is equal 
to 2.37 x 10~ 21 cm 2 /H. For constant A gr (0), the solution of 
eqn. JTDJl is 



no{t) = no(0)e t/ta 



(12) 



where no (0) is the initial gas-phase oxygen atom abundance 
and t acc is their accretion time-scale 



^acc — 



1 /27rm \ 1/2 
A BI (0)n H S \ kT ) 

Integration of eqn. (j4} then yields 
d(t) = d max [l - e -'/*-°] 



(13) 



(14) 



where d ma x = (18mHno(0)/A gr (0)nH/9icc) is the maximum 
thickness of ice that can be attained when all the oxygen has 
frozen out. With d m set from the observational data through 
equation J5), the mantle formation time is found from 



tic 



tacc In 



dn 



(15) 



In order to calculate ti cc in the Coalsack, we adopt a 
sticking coefficient of unity, and assume a gas and dust tem- 
perature of 10 K (Jones et al. 1984). Experimental mea- 
surements of amorphous water ice deposited at low tem- 
peratures indicate densities in the range 0.6-1 g cm -3 , de- 
pending on the details of the ice deposition (e.g. Jenniskens 
& Blake 1994; Stevenson et al. 1999). Therefore, we adopt 
pice = 0.75 g cm -3 . Values of ATjce, Ay, and wh for each of 
the three globules are taken from the observations of Kato 
et al. (1999) and Smith et al. (2002). The only remaining pa- 
rameter required to calculate ii cc is the initial oxygen atom 
abundance in the gas, which depends on the elemental abun- 
dances in the Coalsack, and on the original chemical con- 
ditions in each globule (i.e. how much oxygen was locked 
up in CO and other molecules). For the chemical models 
described in the following section, O atom abundances of 
no(0)/riH = 1.6 x 10 -4 and 2.9 x 10 -4 are assumed, depend- 
ing on whether the chemistry in each globule was initially 
appropriate for translucent (T) or diffuse (D) material. 

Table 1 lists the physical conditions observed in each 
globule, and the corresponding values of d m and U C c derived 
for each object. In all cases the derived mantle thicknesses 
are significantly less than the maximum possible value, d max , 
showing that only a small fraction of the total oxygen is 
frozen out as water, and that each globule must be much sig- 
nificantly younger than the accretion time-scale. The ti ce val- 
ues derived are consistent with these globules having formed 
recently (cf. ~2x 10 5 years; Lada et al. 2004). The princi- 
pal uncertainties in our derived values of ti ce come from un- 
certainties in the values of A gr (0) and the initial gas-phase 
oxygen atom abundance. For example, if we truncate the 
grain size distribution at a_ = 100 A instead of 50 A, this 
reduces the effective grain surface area by a factor of almost 
two-thirds, and thus increases the derived values of ti cc by 
50 per cent. Similarly, because the initial values for no differ 
by a factor of two between models D and T, the derived ages 
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differ by a similar factor. Hence, determining the absolute 
age of each core requires knowledge of the initial chemical 
conditions and the grain size distribution. However, it is 
likely that each globule coalesced from material with simi- 
lar properties, and so the relative ages will be much more 
accurate. 



4 CHEMICAL MODEL 

We now want to use the derived ti cc values to find chemical 
discriminants of the different ages. The values of U ce measure 
the time since the ice mantles began to be laid down, and 
it is well known that mantle formation does not occur until 
the visual extinction reaches a critical threshold of Ay ~ 
3 mag. (e.g. Whittet et al. 1988). In the Coalsack, Smith et 
al. (2002) derive an extinction threshold in the range 2.6- 
7.6 mag. Therefore, the appropriate initial conditions for 
the chemical model at t = are determined by the chemical 
state of the cloud when Ay « 3-8 mag., which are in turn 
determined by the physical evolution of the cloud from the 
diffuse phase to a dense core. In the limit that core formation 
occured rapidly, the initial abundances will be identical to 
those in the diffuse ISM, i.e. principally atomic oxygen and 
nitrogen with carbon and sulphur in ionic form. On the other 
hand, if the cloud evolved slowly, so that it spent a long time 
in the 'translucent' phase (1< Av£ 3 mag.), then the cloud 
material will have been mainly molecular (CO and N2, with 
only O atomic) when Av reached the critical value and ices 
began to form. Therefore, we have considered two models 
corresponding to initally diffuse and translucent chemical 
abundances: the abundances at t — for both models are 
listed in table [5] 

The gas-phase elemental abundances are taken from 
Savage & Sembach (1996), except for sulphur which is as- 
sumed to be significantly depleted relative to the diffuse ISM 
(Ruffle et al. 1999). In our 'standard' model runs we also as- 
sume complete depletion of heavy elements, although the 
effects of non-zero metal abundances are discussed below in 
i|5.1l The chemical reaction network is described in Rodgers 
& Charnley (2001), and has been updated to incorporate 
several recently measured dissociative recombination chan- 
nels (e.g. Geppert et al. 2006). All neutral species except H2 
and He freeze out onto grain mantles with a sticking coeffi- 
cient of unity. We adopt a temperature of 10 K and a cosmic 
ray ionization rate of 3 x 10~ 17 s _1 . Values for tih and Ay 
are set equal to those derived for each of the three globules 
(table CO . 



5 RESULTS 

Figure [1] shows the chemical evolution in Globule 1 for both 
the initially diffuse and initially translucent model runs. Also 
indicated on the graphs are the mantle formation timescales, 
ti cc , derived from the H2O ice observations. The values of t 1CB 
are different for each model since U C c depends on the initial 
abundance of atomic O. Figure [2] shows the evolution in the 
denser Globule 2. Globule 3 has a density intermediate to 
Globules 1 and 2, and the chemical evolution in this source 
is similar to that in Globule 1. 

Comparing figures [T] and [2] it is apparent that although 



Table 2. Initial abundances relative to H2 



Species Abundance 





Model D 


Model T 


c+ 


2.8(-4) 


2(-8) 


c 





2(-5) 


CO 





2.6(-4) 





5.8(-4) 


3.2(-4) 


N 


1.6(-4) 





Ni> 





8(-5) 


s+ 


l(-7) 


l(-7) 


He- 


0.2 


0.2 



Model D corresponds to the abundances expected if the glob- 
ules condense rapidly from diffuse cloud material. In model T 
the initial conditions represent those in low-density translucent 
(Ay Kb 2) gas. a(— b) represents a X 10 _b . 

the chemical evolution of both globules is similar in many 
respects, the difference in density results in a slightly differ- 
ent chemistry in each source. If the globules formed rapidly 
from diffuse cloud material, it takes a long time for the con- 
version of the initially atomic material into predominantly 
molecular form, but in the denser Globule 2 the CO/C ratio 
becomes > 1 after 3 x 10 5 yr, whereas in the lower density 
Globule 1 the transition takes almost 1 Myr. The trans- 
formation of atomic nitrogen into N2 takes even longer in 
each source. The large abundances of C° and N° act to sup- 
press the abundances of species such as CN, NO, and SO, 
which in turn tend to reduce the abundances of many other 
species below their 'typical' dark cloud values. In both glob- 
ules HCN, NH 3 , HC3N, and CH3CN remain at low abun- 
dances for t < 1 Myr, although their abundances become 
somewhat larger in Globule 2 at later times. Based on our 
values of ti ca , Globule 1 is older than Globule 2, and so the 
difference in abundances between the two regions will not be 
as pronounced as if both were the same age, but we predict 
that in Globule 2 the abundances of H 2 CO, HCN, HC 3 N, 
and CS, should be « 2-3 times larger than in Globule 1. 

The alternative scenario, in which CO and N2 formed 
prior to the cores reaching the critical extinction for water 
ice formation to occur, results in larger abundances of sev- 
eral molecules, especially sulphur-bearing species. For exam- 
ple, SO, SO2, and H2S are much more abundant in panels (b) 
of figures Q] and [5] Again, chemical differences between the 
two sources arise primarily due to their different ages, and 
we expect that carbon chain species such as C3H and C4H 
should be 2-3 times more abundant in the younger Globule 
2, whereas NH3, SO and SO2 should be more abundant in 
Globule 1. 

In terms of discriminating between the two model sce- 
narios, there are large differences in the abundances of sev- 
eral species. If the globules formed directly from diffuse ma- 
terial, the abundances of carbon-chain molecules such as 
C4H will be larger than if the gas was initially molecular. 
On the other hand, the reduced C° abundance in the lat- 
ter case leads to larger abundances of many other molecules 
such as H 2 CO, NH3, HCN, HNC, and especially SO and 
SO2. Therefore, observations of these species in the Coal- 
sack can in principle be used to determine which model is 
correct. In order to allow comparison with observations, ta- 
ble [3] lists the column densities for selected species in each 
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Table 1. Globule properties and ice mantle formation time-scales in the Coalsack 



Globule 


N L 


A Y 


"h 


x(CO)* 


dm 


*Ice(D) 


*ioe(T) 




(10 17 cm" 2 ) 


(mag.) 


(10 3 cm" 3 ) 


(io- 5 ) 


(A) 


(10 5 yr) 


(10 5 yr) 


1 


9.1 


11.2 


4.0 


7.5 


85 


2.0 


3.7 


2 


5.7 


13.3 


8.2 


13.7 


45 


0.5 


1.0 


3 


2.0 


9.9 


5.4 


8.6 


21 


0.4 


0.7 



t Smith et al. (2002). i Kato et al. (1999); for globule 3 we have used values quoted for their core 4. ti cc (D) and tj ce (T) are the ice 
formation time-scales corresponding to different values of the initial O atom abundance (see text). 



Table 3. Calculated molecular column densities in each of the three globules. 



Molecule Globule 1 Globule 2 Globule 3 





D 


T 


D 


T 


D 


T 


fc™(10 5 vr) 
• / iccv j - w y 1 / 


2.0 


3.7 


0.5 


1.0 


0. 1 


0.7 


c 


4.6(17) 


4.1 


15) 


7.0(17) 


6.1(14) 


1.3(18) 


1.0(15) 


c+ 


1.0(15) 


1.9 


14) 


1.9(14) 


9.2(13) 


4.0(15) 


1.4(14) 


CO 


1.3(18) 


7.1 


17) 


1.6(18) 


1.5(18) 


8.3(17) 


1.6(18) 


CH 4 


7.6(15) 


6.8 


14) 


5.3(16) 


7.0(14) 


6.1(14) 


7.8(14) 


H 2 CO 


4.0(15) 


4.5 


14) 


1.4(16) 


1.8(14) 


6.1(14) 


2.0(14) 


H 2 


7.4(15) 


1.1 


16) 


3.5(16) 


1.9(16) 


2.7(15) 


2.1(16) 


OH 


6.6(14) 


5.4 


15) 


5.3(14) 


4.0(15) 


2.1(14) 


3.2(15) 


o 2 


4.6(14) 


6.2 


16) 


6.2(14) 


1.7(17) 


1.5(14) 


1.1(17) 


NH 3 


7.8(14) 


4.1 


15) 


9.0(14) 


3.0(15) 


6.8(13) 


1.5(15) 


NO 


1.6(15) 


4.5 


16) 


1.0(15) 


3.3(16) 


3.1(14) 


2.1(16) 


HCN 


6.4(14) 


1.6 


15) 


4.6(15) 


4.7(14) 


1.0(14) 


2.6(14) 


HNC 


2.9(14) 


7.8 


14) 


1.2(15) 


2.7(14) 


3.4(13) 


1.6(14) 


CN 


3.6(14) 


5.8 


14) 


2.7(14) 


1.2(14) 


8.9(13) 


1.2(14) 


CS 


3.8(14) 


8.8 


13) 


4.8(14) 


7.4(13) 


2.4(14) 


1.3(14) 


SO 


6.7(11) 


9.9 


13) 


2.4(11) 


2.5(14) 


3.3(11) 


1.5(14) 


so 2 


2.2(10) 


6.0 


13) 


3.7(10) 


1.8(14) 


5.3(9) 


5.6(13) 


H 2 S 


1.0(11) 


4.5 


11) 


1.8(11) 


6.0(11) 


1.1(10) 
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5.4(13) 


2.1(13) 
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13) 


3.4(15) 


9.0(12) 


2.6(13) 


2.2(13) 


HCO+ 


4.0(13) 


9.8 


13) 


7.2(13) 


1.4(14) 


3.5(12) 


1.2(14) 


O 


9.7(17) 


3.9 


17) 


2.3(18) 
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CH 
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2.6(12) 


2.8(14) 


3.0(12) 


CH 2 


2.9(13) 


6.7 


ID 


5.3(13) 


2.9(11) 


5.7(13) 


3.1(11) 


C 2 H 2 


1.5(14) 


5.7 


13) 


1.7(14) 


5.2(13) 


3.3(13) 


1.1(14) 


N 2 


4.8(17) 


6.4 
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4.2(17) 
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1.3(17) 


6.0(17) 


NH 


2.1(14) 
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2.7(13) 


2.8(14) 
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CH 3 CN 


6.2(13) 


1.3 


13) 


1.4(15) 


1.6(12) 


1.7(12) 


1.5(12) 


CH 3 OH 


9.6(9) 


4.6 


10) 


1.7(11) 


:.0(10) 


1.1(8) 


1.5(11) 


CH 3 OCH 3 


6.3(6) 


4.. 


2(6) 


3.6(8) 


6.0(6) 


1.2(4) 


8.3(6) 


OCS 


1.3(12) 


4.1 


11) 


2.4(12) 


1.7(12) 
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1.9(12) 


N 2 H+ 


9.7(12) 


5.4 


13) 


7.0(12) 


3.0(13) 


3.3(11) 
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a(b) represents a X 10 b cm -2 . 

Column densities are calculated by multiplying the abundance calculated by our chemical model at time t; CG by the H 2 column 
densities determined from the visual extinction measurements of Smith et al. (2002). 

D and T refer to models with the initial conditions appropriate for diffuse or translucent material respectively (see table 
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Figure 1. Chemical evolution in Coalsack Globule 1. (a) The left panels show the model results for initially diffuse cloud abundances, 
and (b) the right hand panels show the evolution from initially transclucent cloud conditions. The gray vertical dashed lines mark the 
ice mantle formation time-scale, t; cc . The time-scale is longer in (b) because the initial abundance of atomic oxygen is smaller. 
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Figure 2. As figure 1, but for Globule 2. 
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Figure 3. Effects of metals on the chemical evolution in Coalsack Globule 2. As in figures 1 and 2, the left and right panels show the 
results for models D and T respectively, and the gray vertical dashed lines mark the ice mantle formation time-scale, fi cc . Solid lines 
show the results for models with zero metals, and dashed lines show the results with a silicon abundance of 10 — 7 . 



of the three globules at the appropriate time, ti ce , calculated 
for each model. 

5.1 Effects of model parameters 

We have also investigated the effects of altering various pa- 
rameters of our models. As discussed above, we used the el- 
emental abundances of Savage & Sembach (1996). However, 
if we consider the only direct measurement taken toward the 
Coalsack, Jensen, Rachford & Snow (2005) reported a low 
elemental O/H ratio of only 6.5 x 10 -5 toward HD 110432, 
a factor of four less than used in our model (although the 
error bars on this measurement are also consistent with a 
'normal' O/H ratio.) We therefore ran our model with the 
abundances of all elements depleted by a factor of four. In 
the case of initially translucent conditions (model T) the ef- 
fects are relatively minor, but for model D we find that the 
reduced elemental abundances act to speed up the transition 
from the diffuse-cloud-type chemistry to the dark-cloud-type 
chemistry, resulting in enhanced abundances of species such 
as H 2 CO, NH 3 , HCN, SO, and S0 2 . However, a key as- 
pect of this change is that it also increases the values of 
tice by a factor of four, since the O atom freeze-out rate is 
reduced. For the observed ice thicknesses in Globule 1 this 
implies that that this source must have an age comparable 
to or greater than the accretion time-scale, which is equal to 

1.2 Myr. This conflicts with other evidence that the South- 
ern Coalsack is a relatively young region (see § 2), and so we 
consider it likely that the elemental abundances (or at least 
the oxygen abundance) in the Coalsack are close to standard 
interstellar values. 

We also investigated the effects of a non-zero metal 
abundance, by running each model with an initial abun- 
dance of Si + equal to 10~ 7 . We find that the effects on 
the chemistry of Globule 1 are minor, but in the denser 
Globule 2 the effects are more important. Figure [3] com- 
pares the abundances of those species most affected by the 



metal abundance in models D and T for Globule 2. We see 
that the effects are largest in model D at late times; at the 
early times implied by the observed water ice column densi- 
ties the presence or absence of metals has little influence on 
the predicted abundances. One exception is HC3N, where 
the early-time abundance is depressed in model T when Si + 
is included. 



6 DISCUSSION 

Most models of interstellar dark cloud chemistry assume ini- 
tial conditions appropriate for diffuse clouds, i.e. they as- 
sume that the cloud formation timescale is much less than 
the chemical timescale. For typical dark cloud densities of 
n(H-2) = 10 4 cm~ 3 , the exact choice of the initial condi- 
tions is not so important, as the chemistry rapidly evolves 
from the diffuse cloud conditions into 'standard' dense cloud 
chemistry. Our results show that, for lower density gas such 
as the globules in the Coalsack, the transformation time may 
be much longer, of order Myr. In this case, the time spent in 
the translucent phase is a key parameter of the model, and 
the dense globules retain a chemical 'memory' of the prior 
conditions. Hence, in principal, molecular observations of 
the dark globules can be used to constrain the dynamics of 
their origins. 

Rachford et al. (2001) observed CO and H2 toward the 
field star HD 110432, and derived a CO/H2 column density 
ratio of ~ 3 x 10 -7 , suggesting that the molecular fraction 
in the inter-clump medium is low. However, at such low H2 
column densities the overall CO/H2 ratio may be low due 
to CO dissociation in the outer layers where Ay< 1, and 
is still consisitent with almost all C in CO in the central 
regions (e.g. van Dishoeck & Black 1989). It is difficult to 
use such observations to tightly constrain the C/CO ratio 
in the interclump medium. 

Although there have been several studies of dust and 
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Figure 4. Integrated pencil-beam H2 column density (i.e. A v \ 
solid line) and water ice column densites (dashed lines) as a func- 
tion of impact parameter for a Bonnor-Ebert sphere with central 
density njj = 10 4 cm" 3 , appropriate for Globule 2 in the Coal- 
sack. Ice column densities are shown at times of 2, 4, 6, 8, & 10 
X 10 5 yr. 



extinction in the Coalsack, there is a dearth of molecular 
line observations of this region. As far as we know, the only 
molecules observed in the Coalsack are CO (Huggins et al. 
1977; Codina et al. 1984; Kato et al. 1999), OH, and H 2 CO 
(Brooks et al. 1976). The latter observed the 6 cm formalde- 
hyde ln-lio line in absorption, and derived a H2CO/H2 
abundance of ~ 10~ 9 , at least an order of magnitude smaller 
than the abundances we derive for the globules in Table 3. 
However, the 6 cm line absorption feature appears to arise 
preferentially from the lower density regions surrounding the 
dense core rather than from the dense core itself (Snell 1981). 
This would appear to be confirmed by the fact that Brooks 
et al. report the emission from the H2CO is extended over 
a size of 12 arcminutes (similar to that of the OH emis- 
sion) , and also by the fact that our chemical model predicts 
a H2CO abundance of ~ 10~ 9 in low-density, translucent 
material. A second effect to consider is that extinction maps 
reveal the densest regions of the globules are likely of order 
1-3 arcminutes in size (Bok 1977; Jones et al. 1980; Lada et 
al. 2004), smaller than the 4 arcminute beamsize of Brooks 
et al. Hence, H2CO features arising from the centers of the 
globules will suffer from beam dilution. We conclude that 
the 6 cm observations of Brooks et al. most likely trace the 
interclump medium. 

We have only considered the case of a single sightline 
through a constant density medium since that is the na- 
ture of the observational constraints currently available to 
us (Smith et al. 2002). In reality, such prestellar globules 
have spatial density gradients, and the observed ice column 
density depends on the particular line-of-sight path integral 
through the cloud. Recent studies of the infrared extinction 
of field stars which lay behind these globules have allowed 
detailed maps to be made of the visual extinction through 
several globules (Alves, Lada & Lada 2001; Kandori et al. 
2005). For isolated globules, the dust density profile is usu- 
ally well matched by a Bonnor-Ebert sphere (Alves et al. 
2001; Racca et al. 2002) and many appear to be on the crit- 



ical limit of gravitational instability against runaway col- 
lapse (Kandori et al. 2005). The ice formation time-scale 
derived for any sightline will always lag that of the dynam- 
ical time-scale over which the mass distribution was assem- 
bled. Hence, if 3 micron observations towards fields stars 
behind such globules can be made to allow their radial ice 
profile to be constructed, this could eventually be employed 
as a good measure of the dynamical time. The analysis pre- 
sented in this paper allows us to calculate the water ice col- 
umn density as a function of both time and offset, £, from 
the core center in Coalsack Globule 2, where nu(r, £) is from 
a Bonnor-Ebert sphere (Racca et al. 2002). As an example, 
Figure (|4]) shows that the ice profile has a distinctive shape 
at each time, and so fitting the observed iVicc(£) distribu- 
tion with such models will allow a global derivation of tice 
for each source. 



7 CONCLUSIONS 

Oxygen atom accretion and reaction to form H2O ice is a 
'one-way' process (at least until a protostar has formed), 
is essentially unaffected by gas-phase chemistry, and begins 
only when the dense core forms. Hence, ice column densities 
can be employed as an independent measure of physical age. 
A combination of ice absorption and molecular emission ob- 
servations can thus lead to quantitative measures of globule 
age, and allow correlations between physical and chemical 
evolution to be explored. 

We have determined the 'ice mantle age', ii cc , of three 
dense globules in the Southern Coalsack, based on the time 
needed to deposit the H2O ice abundances observed by 
Smith et al. (2002) in these sources. In conjunction with 
a simple chemical model, we have used these values of ti ce 
to predict the abundances of numerous gas-phase molecules 
toward each of these sources. These predictions can be easily 
tested by future molecular line observations of the Coalsack. 
Our main findings can be summarized as follows: 

• In each globule only a small fraction of the total avail- 
able gas-phase oxygen is frozen out as water, indicating that 
each must be significantly younger than the accretion time- 
scale. 

• We derive values for ti CB of a few times 10 yr. for Glob- 
ule 1 and < 10 5 yr. for Globules 2 and 3. These ages are 
consistent with the dynamical age for Globule 2 derived by 
Lada et al. (2004). 

• The principal uncertainty in determining exact values 
for iico arises from our lack of knowledge of the gas-phase 
O atom abundance in each source and the grain size distri- 
bution. However, assuming each globule formed from chem- 
ically similar gas, the relative ages of the globules are easily 
constrained. 

• The chemical evolution is very sensitive to the initial 
conditions when the visual extinction becomes large enough 
for water ice to accrete. In particular, whether carbon and 
nitrogen are atomic or locked up in CO and N2 has a pro- 
found influence on the chemistry: if C/CO > 1 at t = 0, the 
abundances of many species are less than if the material is 
initially in molecular form. 

• As the globules age, the abundances of some species 
change, and should be different in each of the three globules. 
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Observations of several species in each of the three globules 
can be used to test the validity of our model. 

• The relative abundances of sulphur-bearing species, 
particularly SO and H 2 S, appear to be useful diagnostics 
of the differences between the various models. 

• Combining water ice observations with gas-phase 
molecular line observations toward pre-stellar cores provides 
a method of dating these objects. Such measurements can 
potentially discriminate between 'fast' and 'slow' scenarios 
for molecular cloud formation. 
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